      SUBROUTINE SETUP
      RETURN
      END
C   HS61
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS61'/
      DATA X/NX*0.D0/
      DATA N/3/ , NH/2/ , NG/0/
      DATA TAU/.1D0/ , TAU0/1.D2/ , DEL0/1.D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,2)/-1,0,0,-1,0,0,-1,0,0/
      END
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=4.D0*X(1)**2+2.D0*X(2)**2 +2.D0*X(3)**2 -33.D0*X(1)+
     1   16.D0*X(2)-24.D0*X(3)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)= 8. D0*X(1) -33.D0
      GRADF(2)= 4. D0* X(2)+16.D0
      GRADF(3)= 4. D0* X(3)-24.D0
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      IF( I .EQ. 2)    GOTO 100
      HXI=      3.D0*X(1)-2.D0*X(2)**2 - 7.D0
      RETURN
 100  CONTINUE
      HXI=      4.D0*X(1)-X(3)**2 -11.D0
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      IF( I .EQ. 2)        GOTO 100
      GRADHI(1)=3.D0
      GRADHI(2)=-4.D0*X(2)
      GRADHI(3)=0.D0
      RETURN
  100 CONTINUE
      GRADHI(1)=4.D0
      GRADHI(2)=0.D0
      GRADHI(3)=-2.D0*X(3)
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
      RETURN
      END

